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Abstract: In this work we study the diamagnetic properties of a perfect 
quantum gas in the presence of a constant magnetic field of intensity B. We 
\ investigate the Gibbs semigroup associated to the one particle operator at fi- 

■ nite volume, and study its Taylor series with respect to the field parameter 

oj := eB/c in different topologies. This allows us to prove the existence of the 
thermodynamic limit for the pressure and for all its derivatives with respect to 



o 

(N 



u (the so-called generalized susceptibilities). 
MSC 2000: 82B10, 82B21, 81V99 

Keywords: semigroup, magnetic field, thermodynamic limit. 

Contents 

m ; . , 

t— I , 11 Intr oduction and results.1 1 
■ ll"l" *l^cIaTion"^ith"tlrc*'de Haas-van Alphcn fdHvA) effect I 4 

12 Ana lyticity of Gi bbs semigroups! 4 

12.1 -Analvticitvl ~~~ 4 

2.2 Proof of Theorem Ol 8 

'"5 ' 12.3 Analvticitv of the semigroup's integral kcrncll 9 

1 ' 

<~| ■ 13 Regularized expansion! 11 

14 Large volume behavior! 17 

\ 15 Thermodynamic limit for magnetic susceptibilities! 22 

JTj ■ 5.1 The proof of Theorem O 26 

rS 

S3. 

1 Introduction and results. 



This paper is motivated by the study of the diamagnetic properties of a perfect 
quantum gas interacting with a constant magnetic field, B := Be$,B > 0, 
e3 := (0,0, 1). The system obeys either the Bose or the Fermi statistics. Since 
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we are only studying orbital diamagnetic effects, we consider a gas of spinlcss 
and charged particles. 

We are mainly interesting in the bulk response i.e. the thermodynamic limit 
of the pressure and its derivatives w.r.t. cyclotron frequency u := eB/c. As in 
|BCL| we use the term generalized susceptibilities to designate such quantities. 

This question has been already addressed by several authors. Thus one finds 
results concerning the existence of the large volume limit of the pressure for both 
Fermi and Bose gases |ACIIABN2| . the magnetization for a Bose gas O MMP 
and the magnetic susceptibility for a Fermi gas [ABN2 . In |BCL| . extensions 
of these results to the case of generalized susceptibilities were announced. 

This paper is the first in a series of two devoted to the rigorous proof of the 
results announced in |B(JL| . Here we consider the regime in which the inverse 
temperature (3 := l/(kT) is positive and finite and the fugacity z = belongs 
to the unit complex disk. Such conditions were also used in |ACIIABN 2 MMP . 
But here we also allow any positive value of the cyclotron frequency lo :— e/cB. 
In a forthcoming paper we will extend these results to some larger z— complex 
domains (in fact to D £ defined below). One can also find different aspects of 
this problem in |MMPI IHoIiol |He3]1 . 

The main part of this work is concerned with a new approach to the mag- 
netic perturbation theory for a semigroup generated by a magnetic Schrodinger 
operator. It extends the results given in |ABN1I IC] and heavily relies on the 
use of the magnetic phase factor. This allows us to have a good control on the 
magnetic perturbation w.r.t. the size of the volume |BCI ICN] in which the gas 
is confined. 

Let us now describe our results. Let A be an open, bounded and connected 
subset of R 3 containing the origin of R 3 and with smooth boundary dA. Set 

A L :={xeK 3 , x/Le A};L> 1. (1.1) 

Here we use the transverse gauge i.e. the magnetic potential is defined as Ba := 
(B/2)es Ax. The one particle Hamiltonian 

H L (w) = ±(-iV -coa) 2 , (1.2) 

is first defined in the form sense on Hq(Al), and then one considers its Friedrichs 
extension RSIfl IRSIV) . Thus we work with Dirichlet boundary conditions 
(DBC). 

It is well-known that H[ j (lo) 1 lo <E K, generates a Gibbs semigroup 

{W L ((3,u;) = e- f3H ^ : > 0} (1.3) 

i.e. for all (3 > 0, Wl((3, lo) € -Bi(L 2 (R 3 ), the set of trace class operators on Hl 
[ACllZj. 

Then for (3 > 0, lo £ R, the grand canonical pressure of a quantum gas at 
finite volume is defined as f |Hu1 IXO lATM^ I 

P L {P, lo, z, e) = ■ Tr {ln(l + ezW L (0, lo)} , (1.4) 

P\ a l\ 

where e = — 1 (e = 1) for Bose (Fermi) statistics. Since lo/2 = inf u{H 00 (lo)), 
then the pressure is an analytic function w.r.t. z on the complex domain D e 
with 

D +1 :=C\ (-oo, -e""/ 2 ], D_! := C \ [e^/ 2 , oo). 
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Let n > 1 and define the susceptibility of order n at finite volume by: 

X { [ L) (fi,uj,z,e) := ^-^(/?,w,2,e); (1.5) 

If n = we set xi (/?,^,2,e) : = PL(/3,oj,z,e). 

Our first result describes the properties of the above defined quantities at 
finite volume, and it is given by the following theorem: 



Theorem 1.1. Let (3 > 0. Then the map M 3 u -> W L (fl, u) € B^L 2 
real analytic and admits an entire extension. For each open and bounded set K, 
which obeys K, C D e , e = — f , +1, there exists an open neighborhood M of the 
real axis such that the pressure at finite volume Pl(/3, u>, z, e) is analytic w.r.t. 
(to, z) on M x /C. Let ui £ K, and |z| < 1. XTien /or n > we have (see (|f .3(1 ): 

We now discuss the limit L = oo. First we define the candidates, Xoo j f° r 
these limits. Recall that the one particle operator Hoo(ui) = \{— iV — wa) 2 on 
i 2 (R 3 ), w 6 i, is positive and essentially selfadjoint on C§°(R 3 ). Denote by 
W oc> (f3, w), /3 > the semigroup generated by H^w). Then W 00 (/3,o;) has an 
explicit integral kernel satisfying (see Section 3): 

Note that the right hand side is independent of x. Let [3 > 0, u> > and |z| < f . 
In view of (|f ,6|l . define: 

Poo 09, w, z, e) := ~ V ^^-G^O, 0; fc/?, u), (1.8) 

which is well defined because of the estimate sinh(i) > t if t > 0. Then by the 
results of jACI IABN2] , we know that 

lim P L (f3,w,z,e) = Poo(P,(jJ,z,e). 

L — >oo 

It is quite natural to choose Xoo^ := d g^n provided that this last quantity 
exists. Note that it is not very easy to see this just from l|1.8|) and j|1.7p . 
We then prove the following 

Theorem 1.2. Let (3 > 0, co > and \z\ < 1. Fix n > 1 and define 

pin p 

X t\l3,uj,z,e):= ^^(13^^,6). (1.9) 
TTien we have the equality : 

x£>(p, = |£ ( T )fca aS° (0 ' 0; fc/3 ' " } - (L10) 

" fc>0 w 



Moreover, 

(i " ) (/?, W ,z,e) = x^ ) 
uniformly on [/?o,/?i] x [wo,wi], < /?o < /?i < oo cmd < < Wi < oo 



lim X r03,w,«,e)=x£ , ()9,w,*,e) (1.11) 

1, — s-oo 
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1.1 Relation with the de Haas- van Alphen (dHvA) effect 



Our results can be easily extended to the case of more general Bloch electrons, 
that is when one has a background smooth and periodic electric potential V. 
More precisely, let us assume that V £ C°°(R 3 ), V > 0, and if T is a periodic 
lattice in R 3 then V(-) = V(- + 7) for all 7 e T. Denote by the elementary 
cell of r. In this case, the grandcanonical pressure at the thermodynamic limit 
will be given by (we work with fcrmions thus e = 1) 

P 00 (/3,cj,z) = ^J2^-T- \k\ I Goo(x,x;fc/3,c)dx, (1.12) 
P fe >-L I I Jn 

where Goo(x, x'; k(3, to) is the smooth integral kernel of the semigroup generated 
by |(— zV + uja) 2 + V. This formula only holds for \z\ < 1, but it can be 
analytically continued to C \ (— 00, — 1], see |AC| or |HeSj| . 

Now one can start looking at the behavior of -Poo (A z) as function of u, 
in particular around the point loq — 0. Working in canonical conditions, that 
is when z is a function of /3, lo and the fixed particle density p, then one is 
interesting in the object 

Poo(/3,w,p) := Poo{(3,uJ,z(l3,uj,p)). 

A thorough analysis of the lo behavior near 0, involving derivatives with respect 
to lo of the above quantity, has been already given by Helffer and Sjostrand in 

ED- 

Alternatively, one can start from the finite volume quantities, and define a 
Zl(P,lo, p) as the unique solution of the equation p^ := /3z9 z Pl(/3, lo, Zl) = p 
and Pl(0, p) := Poo(P, Zl(/3, lo, p)). Is it still true that at large volumes we 
have for example that 

8Zp L (J3, to, p) ~ SZpM lo, p), n > 1 ? 

The main achievement of our paper is that at least for small densities (which 
fix \z\ < 1) the answer is yes. In a companion paper we will prove that this is 
true for all z € C \ (—00, — 1]. 

We end the introduction by giving the plan of this paper. In Section 2 we 
discuss the analyticity of the Gibbs semigroup with respect to lo in the trace class 
sense. The trace norm estimates we obtain depend on the size of the domain, 
due to the linear growth of the magnetic potential. Using magnetic perturbation 
theory we manage to regularize the trace expansions and to extend these results 
to the infinite volume case in Sections 3 and 4. Finally we prove the existence 
of thermodynamic limits in Section 5. 

2 Analyticity of Gibbs semigroups 
2.1 Bi- Analyticity 

Let Al, L > 1 be domains of R 3 as defined in (|l.ll) . In the following we will 
denote respectively by ||T||i, ||T|| 2 and ||T|j, the trace norm in Bi(L 2 (A L )), the 
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Hilbert-Schmidt norm in £?2(L 2 (Az,)) and the operator norm in B(L 2 (Al)) of 
T. 

In this section we study the ui expansion of Wl(/3,u>). This question has 
been already considered |HPI IXBN2 , Z in connection with the B\ analyticity 
of Wl(/3,oj). Combining their result with our analysis below, this gives the 
following. Define the operators: 

R ltL (P,u) :=a-{X7K+ua)W L {f3,u), (2.1) 

■= \a 2 W L {[3,u). (2.2) 

Both operators R\.l{P,uo) and R,2,L(f3,oj) belong to B(Al) and we have the 
following estimate on their norm: 

Lemma 2.1. For all [3 > 0, to > and L > 1, there exists a positive constant 
C such that: 

\\Kl\\ < ^ and \\R 2 , L \\ < CL 2 . (2.3) 

Proof. Let tp e L 2 {A L ). Since W L (J3, u)L 2 (A L ) C Dom(#z,(w)) |Kj, we have 
after a standard argument (note that the absolute value of the components of 
a are bounded from above by diam(Ai) ■ L): 

\\a-{i\7 x + Lja)W L {l3,Lo)ip\\ 2 < C L 2 {H L {u)W L {[3,u)ip,W L {f3,u)<p) 

<^M 2 , (2.4) 

where the last estimate is given by the spectral theorem. The second bound of 
(|2.3(l is obvious. □ 



Remark 2.2. Due to the diamagnetic inequality (see (|2.3()(l below), we have 
for all [3 > and lu e R: 

\\W L (J3,u)\\i = Tr {W L {fl,u>)) < (2 ^ 3 )3/2 - (2.5) 

Then both operators -Ri,L,i?2,L are trace class, since we can factorize the oper- 
ator R 1>L (J3, w) = J2i,zG9/2, u)W L (J3/2, w). 

For n > 1 define: 

V n (fi) := {0 < t„ < r n _! < ... n </3}c R n . (2.6) 

Let (ii,...,i n ) € {1,2}™. Lemma f2 . II allows us to define the following family of 
bounded operators: 



I n ,L(ii,...,in){P,u):= W L (p-n,u)Ri uL (Ti-Ti t u) (2.7) 

JT> n {0) 

■ Ri 2 ,L(T2 ~ T 3 ,U))...Ri n _ llL (T n -i - T n ,Oj)Ri n! L(r n ,Uj)dT, 



5 



where dr is the n-dimensional Lebesgue measure. These operators are in fact 
trace class, and we will estimate their trace norm later. Let n > 1, (ii, ■■■ 1 i n ) € 
{1, 2} n and Xk be the characteristic function, 



Xtih, -,ik) ■= 



1 if ii + ... + i k =n 
otherwise. 



(2.8) 



Then we have 



Theorem 2.3. Fix (3 > 0. Then the operator-valued function R 9 u> h- > 
Wl{/3,uj) £ Bi admits an entire extension to the whole complex plane. Fix 
ujq > 0. For all u> € C we Ziaue 



W L (P,LU) = J2 



{w - w Q ) n d n W L 
nl duj n 



(/3,Wo), 



(2.9) 



1 d n W L 
n! doj n 



09,u;o) = ^(-l) ft ^ xt(k,...,ik)Ik,L(h,.:,ikm^). (2.10) 
ije{i,2} 



fc=i 



Moreover, there exists a positive constant C independent of n > 1, (3 > and L 
smc/i that: 



1 9™iy L 



n! c*w n 



(ft wo) 



< C r 



/?3/ 2 



1 



[(n-l)/4]! 



(2.11) 



for aZZ w e C, {Wl(/3, ui), /3 > 0} is a Gibbs semigroup with its generator given 
by the closed operator Hl(ui). 

Remark 2.4. This theorem implies that the trace of the semigroup Wl is an 
entire function of ui and by (12.9(1 : 



n=0 



duo 11 



(2.12) 



Proof of Theorem We will use here some results from HP] and QVBN2 , 
which we briefly recall. Let uiq > 0. For ui G C set 8uo := ui — luq- Then the 
operator 

H L (u>) - H l (lu ) = (Su>) a • (iV x + w a) + a 2 (2.13) 

is relatively bounded to Hl([3, u>q) with relative bound zero. Note that for /3 > 
from (|2.13|1 we have in the operator sense on Hl 

Rl(P,w) ■■= (H L {u) - H L (u; ))W L (P,u; ) 

= 6u R lfL (p, lu ) + {5uj) 2 R 2 , l ([3, cj ). (2.14) 

For every compact subset K c C, and due to the estimates l|2.3[) . this operator 
satisfies 

/ dr sup ||-Rz,(r, w)|| < oo. 
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Let < Pi < /3q < oo. Then the series 

oo 

W l ((3,u,ljo) = ^2{-l) n wi n) (fl,u,uo) (2.15) 

n=0 

where W { ° ] (f3 , u , u ) = W L (0,u) O ) and for n > 1 
W { L n) (f3,u,u Q )= / dr^^-r,^)^^)-^^))^"" 1 ^^^) (2.16) 



is uniformly B\- convergent on K x [/3\, (3q\. This result was obtained in ABN2 . 
Since Wjj(/3,u>,u)o) is the uniform limit of a sequence of entire Bi-valued func- 
tions, it follows via the Cauchy integral formula that Wl{/3,uj,liJo) is also Si- 
entire in u>. Moreover, for real u> it coincides with the operator c~^ Hl ^\ (3 > 0. 

What we do here is to identify its nth order derivative with respect to u. 
From H2.16 1) and (|2.14|) a simple induction argument yields the following finite 
rearranging: 

N N n 

^(-l)"^" ) (/3,^c ) = T^ i (/3,^) + E(M"E(- 1 )' £ E Xl(k,-,ik) 

n=0 n=l fc=l ije{l,2} 

• h,L{h, —,ik)(f3,uo) +^jv+i,i(/3,w,wo), (2-17) 
where 

n N+hL (l3,u,u ) (2.18) 

2N N 

= E (^rEC" 1 )* E xS(ii,...,ifc)4i(ii»-.ifc)C8»wo)- 

n=iV+l fc=l f J e{l,2} 

Now differentiation with respect to u) commutes with the limit N — > oo, again 
due to the uniform convergence and the Cauchy integral formula. Hence ()2.9JI 
is proved, since the nth order derivative of YljLo(~ l) J Wi (P, w, w ) at lo = u> 
equals the right hand side of Ij2.1t)|l if N > n. 

In the second part of the proof, we use the methods of |ABN2| in order to 
estimate the i?i-norm of the operators Ik,L(ii, ■■■,ik) as claimed in l|2.11|l . We 
first have: 

\\i k , L {ix,...,ik)iP,uo)h < (2.19) 

dr \W L (Ji ~ Tl,Wo)Ri lt L(n - T2,(jJ )~.Ri h ,L(Tk,Uo)\\l- 



Iv, 



Recall that the Ginibre-Gruber inequality read as [XBN2 , 
k / k \ 



Y[AiT(tt)\\i < [UWAW Tr T(t Q + h + ... + t k ) (2.20) 



;=o \i=o 



where {A;,0 < I < k} are bounded operators and T(t),t > is a Gibbs semi- 
group. Then taking A := W L (^,u ), Ai := L (3=g±k iUo ) if I > 1 
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(we put Tk+\ = 0) and T(t) = Wl{^,u)q). On Dk{(3), we have the estimate 
Poll < 1 < Jj=ti and by the Lemma IO for I > 1 



\\M = \\tk u LC- L ^r ± ,m)\\ < const - iti-fiSL. (2.21) 

2 y/Ti - T/ + i 

Let fk ■ T>k(P) — * M defined as 

/fc(r) := Tra v \r r ' (2 ' 22) 

V(P - n)(Tl - T 2 )...(r fe _i - TfcjTfc 

and it satisfies 

Mr)dr=^^. (2.23) 

Let ii + ... + i k = n. Then from l(2~5|). l(2~2T)j> . (|2~2"T|) and lf2~2^|) . we obtain the 
existence of a numerical constant C, such that for every (3 > 0: 

ll4i(<i.-»i*)G9^)lli < • (2 - 24) 

Thus we have the estimate (see l|2.1U[l ): 

p fe=ii 3 -e{i,2} 1 

But a lot of terms in the above sum are zero, since Xfe(*i) — j *fe) = if fc < [^rp 1 ]. 
Since T is increasing, we can give a rough estimate of the form: 

r /*±i) ~ r / [(n+i)/2]+i \ ~ ^ [(n-l)/4]r 1 ' 



fc=l i,e{l,2} 



^ [(n+l)/2] + l 



□ 



2.2 Proof of Theorem fTTTL 

The analyticity properties of the pressure are now easy to prove once we have 
the Bi-analyticity of the Gibbs semigroup. See ABNl] for details. 

Now let (3 > 0,lu > and \z\ < 1. Since \\zWl(P, < 1, the logarithm in 
the pressure at finite volume can be expanded and then: 

P L (J3,u,z,e) = e/(p\A L \)J2(-^) k /kTrW L (k[3,Lu). (2.27) 

k>l 

Starting from Definition (|1.5() . and using Theorem 12.31 we obtain: 

X ^(P,u J , Z ,e) = e m A L \)^ez) k /k Tr ( g^Mf^ ) . (2.28) 

fc>0 ^ ' 

Note that (|2.11(l insures that the growth in k which comes from the trace of 
the nth derivative of WL(kf3,ui) is not faster than some polynomial, but since 
\z\ < 1 the series in k is convergent. This finishes the proof of the theorem. □ 
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2.3 Analyticity of the semigroup's integral kernel 

In the rest of this paper we will only consider Al = (— L/2, £/2) 3 , L > 1. For 
w6R, Hl(lo) is essentially selfadjoint on 

{cp e C\K L ) U C 2 (A L ), 0| S a l = 0, A4> G L 2 (A L )}. 

Let Gl(x, x';/3, w) be the integral kernel of Wl(/3,w) (see e.g. [AC] ). Stan- 
dard elliptic estimates for the eigcnfunctions of Hl(uj), together with the fact 
that e~^ HL ^°> is trace class imply that Gl(x, x'; j3, lo) is smooth in (x, x') £ 
Al X Al. Moreover, Gl(x, x'; /?, cj) = if either x or x' are on the boundary. 

To prove the next theorem, we need the following result from [2], concerning 
the G 1 regularity up to the boundary of the integral kernel. Let j3 > and let 
Goo(x, x',/3) := Gco(x, x', j3, lo — 0) be the heat kernel on the whole space, i.e. 

Goo(x,x';/?) = ^-L^e-^. (2.29) 

Recall that the diamagnetic estimate reads as |ACj : 

|G L (x,x';/3,w)| <G 00 (x,x';/3), (x, x') e A L x A L , lo>0. (2.30) 
Then we have 

Lemma 2.5. Let (3 > and a; > 0. Tften on A^ x A^ we have: 

|(iV x + wa(x))G i (x,x';/3,a;)| < 4g ^(x, x', 8/3). (2.31) 

VP 

where C = C(f3, lo) = c ■ (1 + /3) 5 (1 + u) 3 and c > 1 is a numerical constant. 

This estimate allows us to define the integral kernels of the operators defined 
in 12.2fl ; more precisely, for (x, x') 6 A^ x A^ we have: 

Ri, L (x,x!;P,u)) := a(x) • (iV x + losl(x))G l {x,x' ; j3,uS), 

R 2lL (x,x';P,u) := 1 a 2 (x) G L (x, x'; /?, lo). (2.32) 

Consider the operator Wl(J3, lo) for complex w, defined by a i?i-convergent 
complex power series in Theorem 12 .31 We will now prove that it has an integral 
kernel analytic in lo: 

Theorem 2.6. Let (3 > and fix loq > 0. 

(i). The operator loq) defined in l|2.10|l has an integral kernel denoted 

by (9"Wl)(x, x'; /3, loq), which is jointly continuous on (x, x') G A^ x A^, and 
obeys the estimate 



1 



d n WL 



duo 11 



(x,x';/?,w ) 



<c n (1+^(1+/)^" n>1 (2 . 33 ) 

^3/2 [n=i]| - v ; 



/or some numerical constant c > 1. 

(ii). For a; S C, £/ie operator Wl{/3,lo) has an integral kernel Gl(x, x'; (3, lo) 
given by: 

G^'-,M = f:^^ (^)(x,x';/3, W0 ), (2.34) 

n=0 ' ^ ' 
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where the above series is uniformly convergent on Al X . Thus Gl is jointly 
continuous on Al X A^ and is an entire function of ui. 

Proof of the Theorem \2.b\ Lemma 12.51 obviously implies for (3 > and u> > 

the estimate: 

|J2i, L (x,x / ;/J ) w)| < J L^G 00 (x,x';8/3). (2.35) 

We also have 

L 2 

|i? 2 , L (x,x';/3,^)| < — G oc (x,x';/3). (2.36) 

In the following, we will often use the uniform estimate with respect to the index 
% = 1,2: 



|i? 4 , L (x,x';i,w)| <r i C 1 -^ii±i^G 00 (x,x';8t) < t < /?, (2.37) 

where d :=d{0,u) = 2\/2C(/?,w). 

Let us start by proving (i). Fix /3 > 0,wo > 0, L > 1, and consider the 
operator ik){(3, u>o), k > 1 defined in (|2.7(l . It admits a continuous 

integral kernel, 2fc,z,(x, x'; /?, w ) := h,L(h, ifc)( x > x'; /?, w ) on A L x A L given 

by 

J fclL (x,x';/? s wb) (2.38) 
= / dr / dyGi(x,yi;/3- ri,w )-Ri 1 ,L(yi,y2;n - T2,w ) 

-Ri fc _ 1 ,i(yfc-i,yn;Tfc-i - r fc ,w )^i )c ,L(yfe,x';Tfc,wo) 



where <iy denotes the Lebesgue measure on R 3fc and Dk(fi) is defined in 12.611 . 
Let i\ + %2 + ... + ife = n. Then by using the Lemma 12.51 the estimate 

|Gi(x,yi;0-Ti,wo)| < 8 3/2 ^^G co (x, y i; 8(/3 - n)), < n < /3, 
and H2.37fl . the following estimate holds on Ai x A^: 

|/ fc , L (x,x';/3,c )| <8 3 / 2 (l + ^)( fc+1 )/ 2 C 1 fc L n / fk(r)drf dy 

JD k (p) Jai 

Go^x, y i; 8(/3 - r 1 ))...G co (y fe , x'; 8r fe ), (2.39) 

where the function is defined in (|2.22|) . Notice that by using the semigroup 
property 

/ dyG 00 (x,yi;ti))...G 00 (yfc,x / ;tfc) < / dyG 00 (x,yi;ti))...G 00 (y fc ,x';tfc) 

= G co (x,x';t 1 + ...+t fc ) (2.40) 

Therefore from Ij2.39|l we get 

|/ Jb , i (x,x';ftwo)| < 8 s / a (l + i 9){* +1 )/ a CjL"G 00 (x J x / ;8)8) / A(r)dr. 

" / C fc (/3) 

(2.41) 
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Then Theorem 12.31 together with l|2.23|l and l|2.41|l show that the operator 

,1 kernel satis 



Qq^t-(/3, luq), n > 1 given by (|2.9H . admits a continuous integral kernel satis- 
fying 



d n Wj 



^(x,x';/3, Wo )| < (8M 1/2 C>!i"G 00 (x,x';8/3)^ £ 



i 4j e{i,2} V 2 ; 



for a new constant C2 = Czifi, u>o) '■= tt 1 ^ 2 (1 + P)Ci((3, ljq). Then by mimicking 
the proof of 12.26J) we get from the last inequality 



d n W L / , 
-^-(x.xjAwq) 



where c is a numerical constant. Since Goo(x,x'; 8/3) < fl67r/3)~ 3 / 2 . (|2.42|) im- 
plies the estimate 1)2.331) and proves (i). Then (ii) follows easily from the pre- 
vious estimate since l/[(n — l)/4]! has a super-exponential decay in n. □ 



3 Regularized expansion 

The bounds obtained in the previous section are not convenient for the proof 
of the existence of the thermodynamic limit of the magnetic susceptibilities. In 
particular the bound on d g ^i L (x, x', (3, ujq) given by (|2.42|) is of order L n . Then 
this gives a bound on its trace of order L 3+Tl , while in view of (|1.6|l we need a 
bound which goes like L 3 . 

In this section, we give an improvement of these estimates. In order to do 
that we need to introduce the magnetic phase <p and the magnetic flux fl defined 
as (here x, y, z 6 and e = (0, 0, 1)): 

</>{x,y) := -e- (yAx) = -0(y,x), (3.1) 

fl(x, y, z) := 0(x, y) + 0(y, z) + cf>(z, x) = i e • {(x - y) A (z - y)}. (3.2) 

Note that /Z is really the magnetic flux through the triangle defined by the three 
vectors, and we have: 

|fl(x,y,z)| < |x-y||y-z|. (3.3) 

For n > 1 and x = yo, yi, y„ some arbitrary vectors in A^, define 

n-l 

Fl„(x,yi, ...,y„) := 0(y„,x) + </>(yfc,yfc+i) (3.4) 

fc=0 

n-1 

= 5Z fl ( x ' yfc ' yfc+i )' if n - 2 

fc=l 

and Fl 1 (x,y 1 ) = 0. 

Notice that due to 1)3. 3|l . we have 

n— 1 fc 

|Fl n (x,yi, ...,y„)| < 'y*- 1 ~ y '" yfe ~ yfe+1 l- ( 3 ' 5 ) 

k=l 1=1 
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Let uj > 0. Consider now the bounded operators given by their integral 
kernels on Al x A^, 

i?i :L (x,x';/3,w) := a(x - x') • (iV x + ua (x)) G L (x, x'; (3, uj). 

# 2 ,l(x,x';/?,u;) := 1 a 2 (x - x')G L (x, x'; /?, w). (3.6) 

Then by the Lemma \'2. 51 a straightforward estimate yields to 

|i? 1 x(x,x';Ac J )| < ^^1g oo (x,x',8/3) < 4CiGoo(x, x', 16/3) (3.7) 

for all (x, x') e Ai x A L . Similarly by H2.30j) we have 

|i? 2 , L (x,x';/3, W )| < ^^G oc (x,x',/3) < ^^(x, x', 2/3). (3.8) 

In the sequel for i = 1 , 2 we will use the estimate on x 

|-Ri lL (x,x / ;/?,w)| <C 3 G 00 (x,x',16/3) (3.9) 

where C 3 = C 3 (/3, uj) := 16Ci(/3, w) and Ci is given in (|2~3T|) . 

Notice that (|3.9Jl provides an uniform bound w.r.t. L and j3 near /3 = on 
the operator kernels. This in contrast with the bound on the norm operator 
of Ri t L',i = 1,2 (see section 2.2, (J2.35JI and l|2.36[l ). Using the Schur-Holmgren 
estimate for the operator norm of an integral operator, 1)3. 9|l eventually implies 

\\Ri,L\\ <C 3 , * = 1,2. (3.10) 
Let x 6 A/,. For fc > 1, m > 0, w>0 and /3 > 0, define the continuous function 



dy 



3 = 1 (i 1 ,...,v)6{l,2}i ^O?) ^ 

(i(Fl j (x,y 1 .....yj)))" 1 Gi(x,y 1 ;/3-r 1 ,a;)^ 1 ,i(y 1 ,y 2 ;r 1 -r 2 ,^) 
m! 

•••-Rij_i,L(yj-i,yj;Tj_i - T j) uj)Ri j!L (y j ,x;Tj ) uj) ) (3.11) 

where in the case of m = we set 0° = 1. 

The main result of this section gives a new expression for the diagonal of 
kernel's nth derivative with respect to uj at finite volume. 

Theorem 3.1. Let (3 > and uj > 0. Then for all x € A^, and for all n > 1, 

one has 

1 r) n G n 

fe=l 

and moreover, uniformly in L > 1; 

\WZ L (x;f3,uj)\<c(m,k) [ +P ^ 3/2 (l + ^-^H 2 , (3.13) 

where c(m, k) — c m+k Yl^-i ^-jr and c is again a numerical factor. 
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Proof. We first need to introduce some new notation. Fix luq > 0. Let lu € C, 
8lu = lu — ujq and Ri^ifi, u), luq), i — 1> 2, Wl{P, u), u)q) be the operators on 
defined via their respective integral kernel given on x A^, by 

^(x.x'^.w.wo) = e 4<5 * x,) i? 4 , L (x,x';/3, Wo ), z = 1,2, 
W L (x,x';/3, W ,w ) = e 45 * x ')G L (x,x';/3,c Jo ) (3.14) 
where is defined in (|3.1|) . We also set 

i? L (x, x'; /?, w, lu ) := 5luR 1iL (x., x'; /?, u, uj ) + (5lu) 2 R 2 .l{* 1 x'; (3, u, w Q ). 

Except a phase factor the kernel of Wl and Ri y L,i = 1, 2 is the same as the one 
of Wl, Ri,L, i — 1,2 respectively. Then they satisfy (|2.30|i . I|3.9() respectively. 
Hence by the same arguments as above, they are bounded operators and 

< 1, H^iH < C 3 (3.15) 

(see $nfy ). Notice also that since \\W L \\ l2 = \\W L \\i 2 < then by fity 

Ri,L,i — 1,2 as well as Ri t L,i = 1,2 are in the Hilbcrt-Schmidt class and for 
(3 > 0, lu > and lu e C, 

£3/2 

ll-Ri.iOMIk. H^ill^OS.w.wb) < Cs||W£(16Aw)||i a < C 3 (27r/3)3/4 . (3.16) 

where C3 = Ca(/3, luq) was first introduced in (|3.9(l . We now define the following 
family of bounded operators on Hl- Let k > 1, {ii, € {l^}*, /3 > 0, 
Wo > 0. For all w G C set 

lfc,z,(ii, ...,ifc)(/3,w,w ) := / (ir Wz,(/3 - ti,cj, cjo)-Rh,£(ti - t 2 , w, w ) 

JD fc (/3) 

■ ■■R lk _ 1: L{Tk-l ~ T k ,UJ,UJo)Ri ht L(Tk,U,U}o), (3-17) 

and for n > 1 

n 

^(A^wo) Xk(iu-,ik)Ik, L (h,-,ik)(0,u>,oj o ) (3.18) 

fe=l ij£{l,2} 

where x]J is defined in 1)2.8(1 . 

Lemma 3.2. Xei N > 1, /3 > 0, u>o > 0. For all lu G C set Sui = lu — ojq- Then 
as bounded operators we can write: 

N 

W L {f3 1 Lu)^W L {(3,Lu,Lu ) + J2 ( - Suj ^ Wn ' L ^^^^ + ( 3 - 19 ) 

n=l 

where TZjy, j l (/3,cj,cjo) a»id l(/^, w , w o) ar e t/ie following bounded opera- 

tors on Hl-' 

2N+2 

n=N+l ij£{l,2} 
• / dr W L ([3 - n, u})U il!L (Ti - T 2 ,LU, LU )...ji iN+uL (T N+1 , LU, LU ) (3.20) 
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where Z?jv+i(/3) is given in (|2.(i|) and 

KP +hL (J3,u,uo) (3.21) 

2JV N 

= E (M w £(-i)* E xUh,..,i k )-i k Ah,.-,ik)(p,u,oj ). 

n=N+l k=l ije{l,2} 

Proof of the lemma. In this proof we fix wo > and omit everywhere the loq 
dependence. We first note that Wl{P, w) is strongly differentiable with respect 
to f3 > (see and satisfies 

By using Proposition 3 from |C] , we can write the following Dyson-type integral 
equation: 

W L {f3,co) = W L (0,uj)- cItW l ((3-t,uj)R l (t,ou). (3.22) 



The above integral is a Riemann integral and converges in the operator norm 
sense. By iterating (|3.22|) we obtain: 

N r/3 r/3-Tt ,/9-tl-...— r»-i 

W L ((3 1 uj) = W L ((3,u;) + y2(-l) n dn dr 2 ... dr n 

Jo Jo Jo 



n=l 



■ W l (P-tl - ... - T n ,u)R L (T n ,w)...R L (n,u) +n$ +1Ji (/3,u) (3.23) 
where 

rP rP—Ti rP—TL—...—Tff 

n$ +1>L (P,Lo) = (-l) N+1 drj dr 2 ... dr N+l 

Jo Jo Jo 

•W L (p-n - ... - t n+1 ,u) ■ R L (T N+1 ,uj)...R L (n,uj). (3.24) 

Then a straightforward change of variables in the integrals of the r.h.s. of the 
last two formulas, yields to 



W L (P,U>) =W L (P,U) + EH-)" / dTW L (P-T 1 ,Cj)R L (T 1 -T2,t 
n=l J D n (fi) 

■■■RhiTn-1 — T n ,LU )R L (r n ,Lu)+n { ^ +hL (f3,u;). 



with 

K < $ +lL (P,w) = (-l) N+1 I cLtW l {I3 - n ,u) ■ R L {Ti,u)...R L (Ttf+i,u) 

Jd n+1 (P) 

(3.25) 

where D n (j3) is defined in l|2.6|l . Recall that Rl = SujRi.l + (Slu) 2 R 2i l- So 
(|3.25(l gives (13.20(1 and a simple induction argument finishes the proof of the 
lemma. □ 
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Continuing the proof of Theorem lff.il From Theorem 12.61 we know that for 
x e Aj,, and /3 > 0, C 3 uj —> Gl(x, x; (3, u>) is an entire function. 

In order to prove l|3.12[l . we will show that for all x e Al we have: 

JV n 

G L (x, x; 0, w) = G L (x, x; (3, u> ) + ]T VF fc "- fc (x; /3, wo) 

«=i fe=i 

+ ^tv+i,l(x;/3,w,w ), (3.26) 

where the remainder term satisfies the property that its first JV derivatives at 
ujq are zero. 

By rewriting Lemma 13.21 in terms of the corresponding integral kernels, and 
looking at the diagonal of these kernels, we have (remember that <^>(x, x) = 0): 

N 

G L (x,x;f3,u) = Gl(x,x;/3,w ) + ^(M" W n , L (x, x; 0, uj, uj ) 

n=l 

+^jv+i,.E,( x > x ; A w > w o) + ^w+i,i( x ' x 5 Z 3 ' w > w o)- (3.27) 

By construction, the two remainders are smooth functions which remain smooth 
even if they are divided by (5ui) N+1 , see formulas H3.20|l and (|3.21(l . This means 
that their first JV derivatives at ujq are all zero. Thus the JVth derivative of 
Gl(x, x; j3, ■) at ujq can only come from the W's. 

What we still have to do is to remove the ui dependence from Ws. For, let 
us show that for 1 < n < JV, x G Al and \Sw\ < 1, 

N 

W niL (x,x;P,u>,u> ) = Y, (ferCifxi^ol + RiVi^^^o), (3.28) 
m=0 

where W™ L were introduced in l|3.11|l . and TZ^ N+1 (x; (3, uj, ujq) has its first JV 

derivatives at ujq equal to 0. Indeed, if we replace the integral kernel of I^l 
from H3.17J1 in the expression of Wn,l from (|3.18(1 . we see that we can add up 
all the magnetic phases, and obtain a factor of the type: 

exp{</>(x,yi) + 0(yi,y 2 ) H V (f>(y k -i,yk) + <A(j/fc,x)}. 

Then this exponent will equal the magnetic flux defined in l|3.4fl . plus an addi- 
tional contribution <J>(x, x) which is zero due to the antisymmetry of the magnetic 
phase. Now if we expand e^ <5w ^ F1 ^ x, " ,yj - ) in Taylor series up to the JVth order, 
we obtain (|3.28|) where the remainder has again the property that its first JV 
derivatives at uiq are zero. Now introduce (|3.28|) in (|3.27|l . and after some alge- 
bra involving the multiplication of two series, we eventually get (|3.26|l . Then we 
can identify the JVth derivative at wo of the kernel's diagonal as the coefficient 
multiplying the JVth power of Suj. The identity l|3.12(l is proved. 

Now let us prove the second part of the theorem, i.e. the estimate Ij3.13|l . 
which is also linked to the natural question "why is formula (|3.12|) better than 
the one from l|2.34|l "? The answer is that W™ L (x; j3, ujq) does not grow with 
L, and we will see in the next section that it even converges when L tends to 
infinity. Let us show here its uniform boundedness in L. 
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Looking at its definition given in (|3.11|) . and using the estimates from (|3.9[) 
together with the diamagnetic inequality, we see that we need to estimate 



|W&(x;/3, W )| <C$J2(-V j E \^ ; > bO / dr f 
3=1 (u,..^)e{i,2 } , ' Jl >- ' J - 

|Fl i (x,y 1 ,...,y i )| 



dy 



Goo(x,yi; 16(/3- Ti))G 00 (yi,y 2 ; 16(n - r 2 )) 

m! 

•■■ G oc (y i _i,y i ; 16(r,-_i - Tj)) G oc (y i , x; 16r,) (3.29) 

Let a = 16m, and identify x = yo. In view of Ij3.5|l and the explicit form 
of the heat kernel, for 1 < I < j — 1, 1 < I' < I; yo,yi, ■ ■■,yj S A^ +1 and 
{ T ij T 2, ■•• r j} G Dj(f3) we need the straightforward estimate 

\yi'-i -Yi'Wyi - y;+i| exp(- J*^ 1 ^ )... exp(- |y 2 fc ~ T *' ) 

< 2a/3exp(-Jf^ 7 )...exp(-^). (3.30) 

Thus (ESIl and (|3~3U|) imply 

|Fl J (x,y 1 ,...,y,)| m G co (x,y 1 ;16(/3-T 1 ))...G co (y„x;16r,) 

i-i 

< (^/a/3r2 3 ^'+ 1 )/ 2 G 0O (x,y 1 ;32(/3-T 1 ))...G co (y„x;32r,) 

< (8j 2 m/3) m 2 3 W+ 1 )/ 2 G 00 (x,y 1 ;32(/3-T 1 ))...G 00 (y,,x;32r,). (3.31) 

Integrating over the spatial coordinates, using the semigroup property 12. 41) [I . 
and then integrating over r variables, one eventually obtains the uniform upper 
bound in L given in Ij3.13|l . 

□. 



Remark 3.3. TheoremWl\ gives us what we need for the purpose of this article. 
One can show that our analysis can be applied in order to get the off- diagonal 
terms of the integral kernel i.e. 9 d ^ n L (x, x'; j3, uj), (x, x') £ x A^ and [3 > 
0,w > 0. In that case we get 



EE(-!) J / dr xU^->ij) I dy 



k=ij=i 



iV 



Dj< - P) (ii,...,ij)e{i,2p 

(»(ECi fl ( x . ym) + fl(x, y 3 , x') + 0(x, x'))) ' 

(n-k)\ 

Gi(x,yi,/3 - T 1 ,aj ) R iuL (y 1 ,y 2 ,T 1 - t 2 ,uj) ...Ri ,,z(yj, x', Tj, i 



(3.32) 
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Remark 3.4. Let j3 > and u>o > and uj E C. From (|3.16[) we have 

\\W L (P - t,uj,uj )R l (t,uj,ujo)\\i < || W L {(3 - t, uj, ujq)\\2\\Rl(t, uj, wq)\\ 2 

C(13,uj,uj )L 3 
~ [fi - r) 3 / 4 r 3 / 4 ' 

where C{(3,uj,ujq) = C^{[3, wo)(| <5w | + | (<W) | 2 ) and C^{P,ujq) is given in (I3.9|l . 
T/ien i/ie Bi-operator valued function re (0, /3) — > (/3 — r, lu)Rl (t, w, Wo) is 
Bi-integrable. Denote by Ul{(3, uj, uJq) — J Wl(P — t,uj)Rl(t,uj,ujq), 

r0 i 16fV 3 

||D' L ^ W , Wo )|| 1 <GL»jf i ( ^_ t)3/4t3/4 <— . (3.33) 

TTie Duhamel-type formula (|3.22|l i/ien implies that Wl(/3,uj) is of trace class as 
a sum of the two trace class operators, Wl and Ul- Consequently the operators 
Ik,i,{iii •••! ik){P, ujo) defined in (|3.17|l are of trace class because the integrals 
only involve B\-integrable functions. 



4 Large volume behavior 

For further applications in Section 4 we need to have a similar result as in 
Theorem 13.11 but with L = oo. The results of Section 2 cannot be applied to 
this situation. On the contrary, we will show in this section that Theorem 13. II 
remains true even if we take L — oo, and the quantities at finite volume converge 
pointwise to the ones defined on the whole space. 

Recall first that the explicit form of the integral kernel of e~^ H °°^;l3 > 
0, uj > is given by 

G^x';/^) = 1 , jffig— e*'**.*') (4.1) 



■ exp 



where the phase 4> is defined in (|3.1(l . 

We start with a technical result. For any x 6 A^, we denote with d(x) := 
dist(x,9A L ). Let M := {(x,x') e A L x A L : d(x) < 1 or d(x') < 1}, and denote 
with xm the characteristic function of M. 

Theorem 4.1. Let (3 > and w > 0. Then for any (x, x') 6 A^ x A^ we /lave: 
|G L (x,x';/3,^)-G 00 (x,x';/3,^)| <2 X m(x,x')G 00 (x,x';/3) (4.2) 
+ C 4 (1 - xm)(x,x')G oc (x,x'; 16/3)e- ( ^i i+ii W I) ) 



| (-iV, - wa(x)) [Gi (x, x'; /?, w) - G^ (x, x'; /?, w)] | (4.3) 

< -%xm(x, x')G 00 (x, x'; 8/3) + C 6 (l - X Af)G co (x, x'; 16p>- ( ^+4^). 
VP 

w/iereC 4 = C 4 (p» = c(l + /?) 6 (1 + w) 4 , C 5 = C 5 (p» = c(l + /3) 5 (1 + w) 3 , 
Cg = Ce(ft, uj) — c(l + /3) 8 (1 + w) 5 <wid c > 1 is a numerical constant. 
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To prove the theorem, we need the following lemma. 

Lemma 4.2. Let (3 > 0, co > and a £ {0, 1}. Then for every (x, x') e x 
we have: 

fl£ G £ (x,x';^w) - fl£ Go, (x, x'; /?, w) (4.4) 

\ ( dr f da(y)d^G 00 (x,y;T,u)[n y -'VyG L (y,y!;p-T,u})], 
1 Jo JdA L 

where do~(y) is the measure on 8Al and n y is the outer normal to 8Al at y. 

Proof. Let f3 > 0, a; > 0. Recall that both Green functions Goo(x, x'; r, w) and 
Gl(x, x'; r, w) satisfy in A^ x A^ in distributional sense the equation 

(i) a T G(x,x';r) = -iHV x -^a(x)] 2 G(x,x';r); r>0 (4.5) 

(ii) G(x,x';t = 0+) = <J( x -x'). (4.6) 
For < r < (3 and on Ai x A^, define the following quantity: 

Q(x,x';/3,r) := f dyG L (y , x; - r, ^ (y, x' ; r, w) . (4.7) 

Then by 14.5(1 it is easy to see that 

a T Q(x,x';/3,r) (4.8) 
= ^jf dy { HV y - ^a(y))2G L (y, x; (3 - r, (y, x'; r, w) 

- G L (y,x;/3-T,w)HV y - wa(y)) 2 G 00 (y, x'; r, w)} . 

Since Gl(x, x';t, w) = if x e <9Al or x' S dAi,, integration by parts gives: 

<9 t Q(x,x';/3,t) = -- / rfo-(y)n y • V y G L (y,x;/3 - T,Lj)G 00 (y,x';T,w). 
1 JdA L 

(4.9) 

Now by integrating with respect to r from 0+ to /3_, and using 1(4.6(1 . we obtain: 
G L (x,x';/3, W )-G 00 (x,x';/3,c i ;) (4.10) 

= \ I dT I d(J (y) n y ■ V y G L (y,x;/3 - r, w)G oc (y, x'; r, ^)- 

Now using the self-adjointness property of the semigroup we obtain G(x, y; r) = 
G(y, x;r), thus we can rewrite 14.10(1 as: 

G L (x,x';/3, W )-G 00 (x,x';/3, W ) (4.11) 
= \ I dT I dcr (y) G oo(x,y;r,w)[n y ■ V y G L (y,x';/3-r,w)]. 

1 JO JdA L 

The lemma now follows from 1(4.11(1 . □ 
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Proof of Theorem \4-l\ Let (3 > 0, uj > and suppose first that (x, x') € M. 
Then l|4.2l) follows from the diamagnetic inequality 1|2.30[) . Let us show (|4.3[1 in 
the same case. We know from (|2.31ll that 

|(-iV x -wa(x))G i (x,x / ;/3,o;)| < -£= G oc (x, x'; 8/3). (4.12) 

VP 

On the other hand, using the observation that -iV x -wa(x) is transformed into 
-iVx-wa(x-x') after commutation with e l( ^ x,x then by direct computation 
from (|4.1|l we get that for all -q > 0, 

|HV x -^a(x))G 00 (x,x';r,/3, W )| < -£i= G °o(x, x'; 2r,/3), (4.13) 

VVP 

where Ci =C[((3,u) = 2(l + r?)(l + w)(l + /3). Then Bl2l and BTHl for r? = 1 
imply H4.3(l . 

Now suppose that (x, x') ^ M. This means that neither points are near the 
boundary. For y 6 <9A then by H2.31|l we have: 

\V y G L fr,x]p-T,u)\ < ■ 7 ^=G QO (y,x,8C9-T)). (4.14) 

VP - r 

By applying the estimates (|4.14l) . i|2.30[l and the Lemma f4. 21 we get 



|G L (x, x'; /?, u) - Goo(x, x'; /?, u)\ (4.15) 

7o JaAi VP - t 

But if y € 9Ai, |x — y| > d(x), then a straightforward estimate shows that for 

_ d 2 (x) 

< f < (3, we have Goo(y, x; 8t) < e ™f> Goo(y, x; 16t). Thus we get 

r.h.s. of (4.16) 

<2 5 d/ dr == / da(y)G 00 (y,x;16(/3-T))G 00 (y,x';16r). 

Jo VP ~T JdK t 

For any t, t' > 0, let us look at the integral 

/ da(y)G co (y,x;t)G 00 (y,x';t')- (4-17) 

JdAr. 



Using the convexity of A^, replacing the integrals on the sides of 8Al by integrals 
on R 2 (thus getting an upper bound) , and using the semigroup property in two 
dimensions, we can show that there exists a numerical constant G > such that 



/ da(y)G 00 (y,^t)G 00 (y,x';t')<C^P^G 0C (^ X ';t + t'). (4.18) 

JdA L Vt VI 

To be more precise, let us look at the integral on the hyperplane defined by 
H : = R 2 + (L/2,0,0): 



/ d^(y)G 00 (y,x.;t)G 00 (y,x';t'), 
Jh 



(4.19) 
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where x and x' are on the same side of R 3 with respect to H . Decompose 
x = xi + X2 and x' = x^ + x 2 where Xi and x^ are the parallel components 
with H, while X2 and x' 2 are the orthogonal components on H . Note that here 
|x 2 | 2 + |x' 2 | 2 > |x 2 - x' 2 | 2 . Since |x - y| 2 = |xi - y| 2 + | x 2 | 2 i f y £ H, we can 
explicitly integrate with respect to y and eventually get Q4.18}! . 
Then we can write: 

/ da(y)G 00 (y,x; 16(/3 - r))G 00 (y,x / ; 16r) < C[ ^ Goo fox'; 16/3). 

Therefore since x, x' satisfy d(x),d(x') > 1 we get 

|G L (x,x';/3,^)-G co (x,x';/3,c i ;)| < 2 i C 1 C' 1 G 00 (x, x'; 16/3) (4.20) 

. g I. 64^ 64,3 > / ,f T IJ 



Due the exponential decay, there are no singularities in this integral, and a 
straightforward estimate gives 14.211 . 

We now use the same method as above to prove (|4.3|l in the case when 
(x, x') ^ M. We know from Lemma T4. 21 that 

(-iV x - Lja(x)) (G L (x, x'; /3, w) - G^x, x'; /3, w)) = ~ / rfr / efo-(y) 

• (-»V X - ^a(a;))G 00 (x, y; r, w) [n y • V y G L (y, x'; f3 - r, w)]. (4.21) 
Then by (I^T^I and B33|l 

|HV x -wa(i))(G 1 (x,x';Aw)-G 0O (x,x';j3,i l ,))j < (4.22) 

4c lC [[ dr[ dff(y) ^(y^-))^(y.;';^). 

Then by using the same arguments leading to (|4.20|l we get 
|(-iV x - wa(i)) (G L (x, x'; /3, w) - G^ (x, x'; /3, w)) | 

< 16 V / ^CiCi 2 G o(x,x';16^)e- ( w f -i«3- ) / dr , , (4.23) 

Jo T (P ~ T ) 

from which l|4.3(l follows. Theorem 14. II is proved. □ 

We now want to prove that the equality 13.12J1 stated in Theorem l3. ll remains 
true even if L tends to infinity. It is well known (see e.g. |ACj ) that for (3 > 0, u > 
and (x,x') e R 3 x R 3 , 

G 00 (x,x';/3,w) = lim G L (x,x';/3,w). (4.24) 

L — >oo 

Our main goal now is to show that this pointwise convergence holds true for all 
the derivatives a G n L , n > 1. 

We need to introduce some notation. Let (3 > and uj > 0. For (x, x') £ 
R 3 x M 3 define 

Ri t00 (x,xf;/3,w) := a(x- x') • (iV x + (x)) Goo(x, x'; /3, w). 
i2 2)0 o(x,x';j9 > w) := i a 2 (x - x')Goo(x, x'; /3, w). (4.25) 
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Le us note that we again have the same type of estimates as in (13.71) , 1)3. 8fl and 
(|3.9I) . i.e. gaussian localization in the difference of the spatial arguments. The 
linear growth of the magnetic potential disappears when one commutes — iV x 
with the magnetic phase, as we have already seen in 1)4.13)1 . 
Now define for x e R 3 , k > 1, m > 0: 

^(x;A.):=^(-l)i ]T \i i; > *i) / drf dy 

i= i (il ,... A .)6{l,2p Jd M Jr3j 

(i^Cx.yx ,-, yj -))) m g 0O (x,y 1 ; < g-r 1 ,^ 1 , 0O (y 1 ,y 2 ;r 1 -r 2) a;) 
to! 

•••-Ri 3 -_i,oo(yi-l,yfe;Tj_i - r j ,uj)R iji00 (y j ,x;T j ,uj). (4.26) 

Since every integrand is bounded by a free heat kernel, and because the flux 
F/j can be bound by differences of its arguments (see (|3.5|) 't . then the above 
multiple integrals are absolutely convergent. Also note the important thing that 
multiplication by |y — y'| m of the free heat kernel only improves the singularity 
in the time variable due to the estimate 

|y - y'|» l e -|y-y'l 2 /* < const • f™/2 e -|y-y'l 2 /(2*)_ (4 27 ) 

The last important remark about W^^x; (3,uj) is that it does not depend on 
x. This can be seen by factorizing all the magnetic phases which enter in the 
various factors of the integrand, and see that they add up to give another Flj, 
which only depends on differences of variables. The remaining factors are also 
just functions of differences of variables. Therefore by changing x we get the 
same value for W™^ after a change of variables (a translation) in all integrals, 
since W™^ only involves integrals defined on the whole space. 
Then we have 

Theorem 4.3. Let (3 > and to > 0. Fix xel 3 and n > 1. Then we have: 
1 r) n f^ 1 f) n C^ n 

^^ ( x,x;A,)=^-^ ( x ) x;A W )=E^(-;A-). (4.28) 

fc=i 

Proof. Fix f3 > and to > 0. Let n > 1 and (x, x') S K 3 . Choose L large 
enough such that x e A^. Then from 1)4.2)) and 14. 3f) we have: 

|G L (x, x';/3, to) - G oc (x, x'; 0,w)\< C^-^er^+^l (4.29) 
|i?i, L (x,x'; P,lj) - i? 1 , 00 (x,x';/3,w)| < C 6 /T V^ 4 ""^, 
|i?2x(x,x';/3,w) - i? 2 , 00 (x,x';/3, W )| < C 4 /T ^e-^+W 1 ). 
Then for all (x, x') g R 3 x E 3 , estimates (I4.29f) show respectively that 
lim Gi(x, x',/3, w) = Goo(x, x',^,w), 

L — >oo 



lim i?j iZ ,(x,x',/3,w) = R ii0o (x,x',(3,uj), 

L — >oo 



for i = 1,2. Then 



lim G L (x, yi ; (3 - n , u;)R h L (yi , y2 ; n - r 2 , to) . ..R i] , L (yj , x; r,- 

L— >oo 
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G OQ (x,yi;/3 - n, w)i2 il>0o (yi, y 2 ; n - r 2 , u})...R ijiL (yj, x; 73, w). 
Furthermore by (|3.9|) and (|2.3U|) we have 

I Gl (x, yi ; /3 - n , , £ (yi , y 2 ; n - r 2 , w) . . .R tj x (y, , x; tj , uj) \ < 

4 3 ^G 00 (x ) y 1 ;16(/3-r 1 ))G 00 (y 1) y2;16(Ti-T 2 ))...G 00 (y J) x;16T j ). 

this last quantity is L-independent and R 3 -integrable by the the semigroup 
property since 

/ dyG oc (x,yi; 16(/3 - n))...G 0o (y i ,x; 16r,) = G 0o (x,yi; 16/3). 

Note that the flux Flj does not influence anything, since it can be bound 
by powers of differences between spatial variables, which will meet the gaussian 
decay of the free heat kernels. Thus they will only affect the time integrals (by 
making them even less singular). 

Then by applying Lebesgue's dominated convergence theorem we get from 
(|3~TT|) and (332) , 

1 f) n C* n 

& "a^ (x ' x ' M = S w t°Z (*; A ( 4 -30) 
00 ■ fc=i 

Now the remaining thing is to show that this also equals ^7 9 g^„°° (x, x, f3, uj). 
Fix w > and choose uj € K such that |<5w| = \uj — uj \ < 1. From the usual 
Taylor formula we can write 

^ f n G 
G l (x,x,0,uj) = ^2(Su) n — -g^(x,x,/3,oj ) 

n=0 

+ ^ )N+1 WTYy. a^ (x ' x '^' Wl) ' (4 ' 31) 

where a->i is between ojo and w. Then by taking L to infinity, we easily get the 
estimate (note that W^ 1 ~ k (x; (3, u>\) is bounded by a constant indepen- 

dent of L, see the estimate from : 

JV n 

Goo (x, x, (3, uj) - Goo (x, x,0,uj o )-J2 (M" E W ^ (*5 ^) 

n=l k=l 

<C(N,P)\5lu\ n+1 . (4.32) 

Since Goo(x, x;/3,w) is smooth in uj (see (EJJl), it follows that the coefficient of 
(5uj) n must equal ^ ^g^ (x, x, /3, w ), and we are done. □ 

5 Thermodynamic limit for magnetic suscepti- 
bilities 

As a consequence of the analysis of the previous section we are now able to 
prove the main technical result of this paper: 
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Theorem 5.1. Let n > 1, P > and lu > 
rfx 



— (x, x;P,u)- (x, x; uj) 



duj r: 



< L 2 C(n,P,u). (5.1) 



where C(n,P,u>) := c(n) — (1 + w) 3n+5 where c(n) on/y depends on n. 

Proof. Let l<J<fc, 1 < & < n, l<m<n — fc, and denote the integrand in 
with: 



_ (iFlj(x,yi,...,yj)) 



(5.2) 



m! 



Goo(x,yi,/3 - Ti,uj)...Ri jt00 (yj,x,Tj,L}). 



Denote also by: 



FTl = ^™ (x, y 1 , • • ■ , yj , Ti , . . . , Tj, p, u) 
(»Flj(x,yi,...,yj)) 



(5.3) 



Gz, (x, y ! , /3 - n , w) . ..Rij }L (yj , x, Tj , w) . 

1 1 hi 

Let n>l,/3>0,w>0, and fix x e R 3 . Then by applying the Theorem 
14.31 we can split the integrals from W's in "inner" and "outer" regions: 



d n G a 
duj n 



-(x,x, J 8,w)=/£(x )J 8,w)+ fl 2(x,/3 > w) 



(5.4) 



where 



TL K> P P 

/r(x,A^):=n!^^(-iy 2 X%ii,-,ij) dr dy 



fc=l j=l 



(u,...» i )e{i,2}J 
i 



( X - yi > •••> yj : n, Tj , p, U), (5.5) 



n k 



5 2(x,/3,o;):=n!^^(-iy £ b)£ / ^ ( 5 - 6 ) 

fc=ij=i (ti,...t J )e{i,2}J i=i jD i(P) 

/ rfyi- / rfy; / dy/+i- / <fyjFjl~ k {x,y 1 ,...,yj,T 1 ,...,Tj,p,u}). 

JR 3 JR 3 \Ar JR 3 JR 3 



Let us now show that 
dx 



A 



f) n G 

f2(x,P,co)-^-^(x,x;P,co) 



j i+/3) 7 "+ 3 



< L 2 f(n,P,u), 



(5.7) 



where f(n,P,uj) := c(n)- — ^ — (1 + u) 3n+2 , c(n) depending only on n. 

From now, for the sake of simplicity we often omit the explicit dependence 
of all variables. In view of (|3.11|l . (|3.12|) an d l|4.26|l . (|4.28|) . we need to estimate 



^ ^ G £ J Ri 1 } L . . .Ri l _ 1 l {Rii ,oo Rii ,L^)Rii-\-i ,00 • --Rij ,00 }• 



(5.8) 



1=1 
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Denote by X (x) the characteristic of {x e A, d(x) < 1}. Thanks to the 
Theorem 14.11 we have 

\(Ri,oc-Ri,L)(x,x-';f3,u)\ (5.9) 
< |a(x - x')(iV x + c^a) {G*, - G L ) (x, x'; 0, u)\ 

<C 7 G 00 (x,x';32/3) ( X (x) + X (x') + e ~( ; w i + ; w i )^ , 

where C 7 = C 7 (/3, w) = c(l + /3) 9 (1 + w) 5 for some numerical constant c > 1. 
But again by the Theorem 14.11 we may use the bound 



|(Goo - G L )(x,x';/3,^)|, \(R 2i00 - R 2 ,l)(x, x'; 0, u)\ < (5.10) 
C 7 G 00 (x,x';32/3) ( X (x) + X (x') + e-C^+T^ . 

On the other hand by (|3.9|l . I|4.13|l and (|4.25|) . the kernel of i£j j0 o, i = 1,2 and 
of i?i,L, i = 1,2 satisfy the inequality 

max{|(iZ i)00 (x,x / ;/?,w)|, |i^ L (x, x'; 0, to)\} < C^G co (x, x'; 16/3), (5.11) 

where C3 = C' 3 0,uj) — c ■ C±, C\ is defined in (|2.37(l and c > 1 is a numeri- 
cal constant which is chosen large enough such that we have Goo(x, x';/3) < 
C£G 00 (x,x';16/?). 

Set y := x. Then (|2~3U|) . JSHH> together with (O and {5~TU)l give 

I T^l I rn 

I^Too " ^1 < C 7 C'f 1 ^rG co (yo,y 1 ;32(/3- n ))... 

G oo (y J ,y ; 32t 3 -) ^(2 X (x ; ) + e"^). (5.12) 

z=o 

Thus from this inequality and 1)3.5(1 , we need to estimate the quantity 

i-i 1 

Q-= (^^|yr-i-yHlyi-y ; +i|) m Goo(yo,yi;32(/3-T 1 ))...G 00 (y 3 ,y ;32T J ). 
1=1 l'=l 

By using l(4.27|l we have 

Q < (8/3j 2 ) m m™2 3 ^+ 1 )/ 2 G oo (y ,yi;64(/3-T 1 )) G 00 (y„ y ; 64r,). 

and then for j < n 

\fZo ~F™l\< 2 3 (" +1 )/ 2 C 7 Cr l (8/3j2) r mm G oo (y , yi ;64(/3- n))- (5.13) 

m! 

Goo to , yo ; 64T,-) £ (2 X ( X/ ) + e - ^ ) . 

By extending the integration with respect to yo, ...y/_i, y;+i..yj on the whole K 3 
space, and using the semigroup property l(2.4()|l and the fact that Goo(x, x;i) = 
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(2irt) 3 / ; 



we get 



A'/ 



- 06763 U + /33/2( n _ fc )! (5 - 14) 



d 2 (* 



<2x(2x(x) + e ) 
for some positive constant c. Moreover, simple estimates show that 
dx(2 X (x) + rw) < C L 2 (1 + y/fj), 



Ar 



where c is also some positive numerical constant. From ()5.5(l and the Theorem 
EH 



dya \ -7^(y ,y ;/3,w) - /£(y ,/3,w) 



(5.15) 



n 



fe=ij=i (ii,...i 3 )e{i,2}j 



/ «fy(i^- k -i^r*)(yo,yi,...,yi,7i,...,r JJ Aw). 



Then (|5.14[) together with l|5.15|l lead to: 



dy [ d J^ (y ,x;/3,o;) - /£(y ,/3,w) 



< L 2 c(n) (1 + ^ 7 " +3 (l + ^) 3 " +2 



where c(n) = (n + l)!c" X)fc=i ^"(n-fc)! ^j— 1 ^~~~~J\ — ~ anc ^ c * s a S am a numerical 
factor. This last estimate clearly implies <|5.7(1 . 

Let us now prove that for all f3 > and ui > 0, g£(yo,/3,w) given in (|5.6() 
satisfies: 



dy gl{y ,fi,uj)\ < L 2 g(n, (3, cj), 



(5.16) 



where g(n,(3,uj) :— c(rt)(l + /3) 7n+2 (l + w) 3 "+ 2 and c(n) is a positive constant 
depending only on n. The same arguments as above leading to the estimate 
(|5"T3|) imply 



l^r°o(yo,---,yj, T i,---,Tj,^,w)| 



(5.17) 



<C 7 C 



<c 7 c 



1-1 IFI,- 



-Goo(yo,yi;32(/3 - ri))...G oo (y 3 ,y ;32T 7 -) 



i-x (8/3j 2 ) m m m 2 3 ^' +1 )/ 2 



Goo (yo , y 1 ; 64 (/? - ri ) ) . . . Goo (y^ , y ; 64r,- ) • 



On the other hand, by the semigroup property (put K C L := R 3 \Ai), 
/ dyi... dy t dy i+1 ... dyj G oo (y , yi; 64 (/3 - n))...G 00 (y J -,yo; 64r.,) 

JR 3 JA= JR 3 JR 3 

= / dy i Goo(yo,y;;64(/3~T i ))Goo(yi,yo;64r i ). (5.18) 
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Then (HSU), (|5~T7|) and imply 



f . n k 

/ dyo<?2(yo,M < m^T^Cf^ 

Ml, , , „._1 



,(5 i+ 3)/2 (8/3j 2 )("- fc )(n-fc)"- fc 



(n-fc)! 



' / / dy ° / d y' G °o(yo, 



y /; 64(/3 - 7i))Goo(yi, y ; 64t z ). (5.19) 



By using the explicit form of the heat kernel given in 1)2.29(1 . a straightforward 
computation shows that 

/ rfyo / rfy;G 00 (y ,y;;64(/3-r ; ))G' 00 (y;,yo;64T ; ) < c^- 

for some positive constant c. Hence we get 



I / dy Q gl(y Q ,p,Lo)\<L 2 c(n)(l+f3) 7n+2 (l + Lo) 3 



7n+2/-i , , ,\3n+2 

A 

where c(n) = n\c n 5Zfc=i J faZfcfezijj — and c is a positive numerical 

factor. This shows 1)5.16)1 . Then 1)5.16)1 and (|5.7|l imply the theorem. □ 



5.1 The proof of Theorem fPl 

We are now ready to prove the thermodynamic limit of generalized susceptibil- 
ities in the grandcanonical ensemble, when the chemical potential is negative 
(fugacity z less than one). 

Let L > 1, P > 0, oj > and \z\ < 1. We know from JTHJ and l(2~2T)> that: 



P L {fi,u,z,e)-PMu,z,e) (5.20) 
e ^ {—ez) k 



/ dx{G L (x,x;fc/3, W )-G co (x,x;fc/3, W )}. 

k>i k 

Then by applying the Theorem 15. II we get 

d n (P L ~ Poo) e x ^ (-ez) k f ( d n G L cTG^\ 

^ = WA g Hr" A/ x I ft^ " "a^ J (x ' x; fcA " } 



In particular, this also shows that the series from (|1.10(l must converge. More- 
over by using again the bound 1)5. l|l in the last formula, we have 



lxi n) -x^l< C W(i + -) 3n+5 ^-E 



\z\ k (l + fc/3) 7 ™+ 8 



Since the series in the r.h.s of this last inequality is finite and L independent, 
this proves 1)1.11(1 . □ 
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